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A diagrammatic theory around atomic limit is proposed for normal state of Anderson Impurity 
Model. The new diagram method is based on the ordinary Wick's theorem for conduction electrons 
and a generalized Wick's theorem for strongly correlated impurity electrons. This last theorem 
coincides with the definition of Kubo cumulants. For the mean value of the evolution operator 
a linked cluster theorem is proved and a Dyson's type equations for one-particle propagators are 
established. The main element of these equations is the correlation function which contains the 
spin, charge and pairing fluctuations of the system. The thermodynamic potential of the system 
is expressed through one-particle renormalized Green's functions and the correlation function. The 
stationary property of the thermodynamic potential is established with respect to the changes of 
correlation function. 

PACS numbers: 71.27.-|-a, 71.10.Fd 



I. INTRODUCTION 

The study of strongly correlated electron systems has 
become in the last decade one of the most active fields 
of condensed matter physics. One of the most important 
models of strongly correlated electrons is the Anderson 
impurity model (AIM)[^1. It is a model of the system of 
free conduction electrons that interact with the system of 
the electrons of the d— or /— shells of impurity atoms. 
The impurity electrons are strongly correlated because 
of strong Coulomb repulsion and they undergo the hy- 
bridization with conduction electrons. This model has 
been largely applied to heavy fermion systems where the 
local impurity orbital is /— orbital ^^1 . 

The interest towards the Anderson impurity model 
has also increased with the advent of Dynamical Mean 
Field Theory (DMFT) within which an infinite dimen- 
sional lattice models can be mapped onto effective impu- 
rity models by using self-consistency conditions''^' ' . This 
model has also been intensively investigated by using the 
method of equation of motion for retarded and advanced 
Green's functions proposed by Bogoliubov and Tiablikov 
f^] and developed in papers t^^^l. 

We propose a diagrammatic method to treat the AIM 
which is an alternative approach to the equations of mo- 
tion and renormalization methods'^'^1 . The first attempt 
to develop the diagrammatic theory for this problem was 
realized in the paper'^"! . These authors used the ex- 
pansion by cumulants for averages of products of Hub- 
bard transfer operators and their algebra. Afterwards, 
other diagrammatic techniques dealing with perturbation 
in the hybridization amplitude while keeping the on-site 
correlation exact have been developed^"'^^"^'^! . Here we use 
the thermodynamic perturbation theory and Matsubara 



[14,15] Qreen's functions, considering the hybridization of 
both groups of electrons as a perturbation. 
The Hamiltonian of the model is written as 
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Ho 



= Ho + H, 
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where Ckcr{C^^) and faifa) ~ annihilation (creation) op- 
erators of conduction and impurity electrons with spin a 
correspondingly, e(k) is the kinetic energy of the conduc- 
tion band state (k, a) , e/ is the local energy of / elec- 
trons, U is the on-site Coulomb repulsion of the impurity 
electrons and N the number of lattice sites. Both energies 
are evaluated with respect to the chemical potential fi of 
the system. The perturbation Hi„t is the hybridization 
interaction between conduction and localized electrons. 
The Coulomb repulsion between impurity electrons is far 
to large to be treated as perturbation and it must be 
included in the main part of Hamiltonian Hq ■ The exis- 
tence of this term invalidates the Wick theorem for local 
electrons. Therefore first of all we formulate the General- 
ized Wick theorem (GWT) for local electrons preserving 
the ordinary Wick theorem for conduction electrons. Our 
(GWT) really is the identity which determines the irre- 
ducible Green's functions or the Kubo cumulants. 

For example the chronological product of four local op- 



erators averaged with respect to zero-order density ma- trix of electrons has the form '^^1 : 



(7/1/2/3/4)0 = (7/1/4)0(7/2/3)0- (7/1/3)0 (7/2/4)0 +5r'''[l,2|3, 4] 



(2) 
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Here the symbol (...)o means thermodynamical average 
on zero order distribution function. The right -hand part 
of equation (2) contains the first two terms which are the 
ordinary Wick contributions and last one named by us 
as irreducible Green's function or Kubo cumulant which 
contains spin, charge and pairing fluctuations of localized 
electrons. In case when the statistical average of oper- 
ators contains 6 operators (7/1/2/3/4/5/g) , ^^^ right- 
hand part contains 3! terms with product of three one- 
particle Green's functions, then there are nine terms in 
the form of product of one-particle Green's function and 
two-particle irreducible function. At last there is a three- 
particle irreducible Green's function ^3 ^*'^[1, 2,3|4, 5, 6]. 

The Green's functions gn [1, ..., n \ n + 1, ..., 2n] which 
appear at the n-th order of the perturbation determine 
the structure of the new diagrammatic technique. Such 
definitions of irreducible Green's functions has been used 
already for discussing the properties of the Hubbard and 
other strongly correlated models [I'^-^i] ^ 

The zero order Hamiltonian of the localized electrons 
can be diagonalized by using the Hubbard transfer oper- 
ators x™" [^^1 . By using these operators it is possible to 
calculate the simplest irreducible Green's functions. 



II. DIAGRAMMATICAL THEORY 

The full one-particle Matsubara Green's function of 
localized electrons in interaction representation has the 
form: 

g„„,{T-T') = - {T,U{r)J^,{T')UiP))l, (3) 

where the index c means connected diagrams. The op- 
erators are taken in the interaction representation , T is 
chronological operator. 

For conduction electrons it is convenient to define the 
local operator 



^H^kCk.- 



(4) 



The corresponding full conduction electron Green's func- 
tion has the form 



C,{T-r') = - {Tb,{T)K,{T')U{P))l, 
where U{(3) is the evolution operator 



rexp(- / Hint(,T)dT) 



(5) 



(6) 



Because the matrix element of hybridization Vk is ab- 
sorbed by local operator ba it is convenient to introduce a 
new parameter A, which will be associated to each vertex 
of the diagrams. In such a way the order of perturbation 
theory will be determined by A and not by the matrix 
element Vk of hybridization which can be present even 
in zero order Green's function. In the last stage of the 
calculation A will be put equal to one. 

In zero order of perturbation theory the Fourier repre- 
sentation of these functions are (ct = —o): 






I- nw 



iuj — €f ibj — Cf — U 
exp(-/3e/) + exp [-/3(2e/ -t- U] 



Zo 
Zo = l + 2cxp(-/3e/)+exp[-/3(2e/ + t/] 



(7) 



Here to = ujn ~ (2n + l)7r//3 are the Matsubara odd 
frequencies. 

For the conduction electrons we have 









|Vk|^ 



k 



E(k)- 



(8) 



The presence in the definition of zero order Green's 
function GjJ of the square of matrix element of hybridiza- 
tion is the consequence of our equation (4) and not of the 
perturbation. 

The thermodynamical perturbation theory gives us the 
results for one-particle Green's functions presented on the 
Fig.l and Fig. 2. The double solid and dashed lines depict 
the renormalized and the thin lines the bar propagators of 
conduction and impurity electrons. The lines connect the 
crosses which depict the impurity states. To crosses are 
attached two arrows one of which is ingoing and one out- 
going. They depict the annihilation and creation of the 
electrons correspondingly. The crosses are the vertices of 
the diagrams and a A multiplier is attached to each of 
them. The index n means {un^Tr^. The summation on 
the index cr„ and the integration on the r„ are intended. 
The rectangles with 2n indices and crosses depict the ir- 
reducible ^i [l,...,n I n -I- l,...,2n] Green's functions. 
The sign of diagrams is determined by the parity (even 
or odd) of the permutation of the Fermi operators neces- 
sary to obtain the diagram. Using Feynman's rules and 
the correspondence above, it is possible to establish the 
next equation for the diagrams shown in Fig.l: 




5r"[l,2|3,4] 



ffril,2,3|4,5,( 



FIG. 1: Diagrams of first orders of perturbation theory for the conduction electron propagator. The thin soUd hues represent 
the conduction electron propagator and dashed lines the impurity electron propagator of zero order. Double line represents the 
full propagator. 
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FIG. 2; The diagrams for impurity electron propagator ga-a-'iT ~ ''"')• The last three diagrams contain the correlation contri- 
butions. Two of them are strong connected and the last is weak connected. 
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G,,,(r-r'|A) - GL'(r - r') + ^ Idn IdrsGL.lr - Ti)Ag.,.,(ri - r2|A)AG:^,,,(T2 - r'). (9) 
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On the basis of the diagrams depicted on the Fig. 2, type equation for g^ 
instead it is possible to establish the following Dyson's 



9aa'{T ~T |A) 



A„ 



'(r - T |A) + ^ dri dT2Kai{T - ri)AG"^^^(Ti - T2)\ga^a'{T2 - T 



(10) 
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where 



Ka'{T-T I A) 



gL'ir 



Z,,-(t-t |A).(11) 



Here Z^^i is the new correlation function which contains 
an infinite sum of the irreducible Green's functions. As 
it was underlined above this function contains all spin, 
charge and pairing fluctuations and is the main element 
of our diagram technique. 

Diagram representation of the correlation function 
^aa'{T — T |A) is depicted on the Fig. 3 

The series expansion of the full propagator Ga-a' iT~T') 
can gives us more detailed representation of this quantity. 

By using the Fourier representation of Matsubara func- 
tions on the base of equations (9)and (10), we have: 



G^iiulX) 



G'A^ 



1 - AAiLo\X)GO{iuj)X^ 



(12) 



The diagrams which determine the mean value of the 
evolution operator ([^(/3))o have not external lines and 
are named vacuum diagrams. Between such diagrams 
there are connected and disconnected ones. The discon- 
nected diagrams can be summed and the result of such 
summation is equal to the exponent of connected dia- 
grams. The result of such summation permits us to for- 
mulate linked cluster theorem. It has the form: 



(C/(/3))o = exp(C/(/3))S, 



(16) 



where (t/(/3))g is the infinite sum of vacuum connected 
diagrams. This quantity is equal to zero when hybridiza- 
tion is absent. 

Therefore thermodynamic potential is equal to 



gaiMX) = 



Aa{iuj\X) 



l-A<,(it^|A)G0(iw)A2' 



(13) 



The equation (12) for conduction electron propagator is 
the Dyson one with mass operator determined by the 
correlation function of impurity electrons: 



^aiiuj\X) = X'^A^{iuj\X). 



(14) 



When A is equal to one these quantities coincide. 

The equation (13) for impurity electrons is of Dyson 
type and coincides, for A = 1, with other equations ob- 
tained for strongly correlated electrons [I'^-^i]^ jj^ gq^a. 
tions (12), (13) the parameter A can be taken equal to 
one and can be omitted. 



III. THERMODYNAMIC POTENTIAL 

The thermodynamic potential of our strongly corre- 
lated system is equal to 

1 



F = Fo- -In ([/(/?)) 0, 



Fo 



llnZo-|^ln[l + exp(-/36(k))] 



(15) 



^ = ^0 - ^ (C/(/3))S , 



(17) 



In Fig. 4 are depicted some of the simplest vacuum 
diagrams. The first three diagrams are of chain type 
and are originated from the ordinary Wick contributions. 
The last three diagrams contain the correlation functions 
and are determined by the new contributions of GWT. 
The factor — , where n is the perturbation theory order, 
present in these diagrams makes it difficult to carry out 
the summation over n. As is usual in such cases '^^l we 
employ a trick, that of integrating over the interacting 
strength A. The result of this procedure is depicted in 
Fig. 5. Now we shall use the diagrams of the conduction 
electron propagator Ga depicted on the Fig.l and these 
of Act from Fig. 3 to combine them in such a way to obtain 
the vacuum diagrams of Fig. 5. Either diagrams of Fig. 5 
of the n order of perturbation theory can be considered as 
the product of the contribution of order ni from Go- and 
of the contribution of order n2 from A^- with the condi- 
tion that ni+n2 = n. There are in general case different 
possibilities to arrange such contribution and the number 
of these possibilities is determined by the numerator of 
the fraction before the diagrams of Fig. 5. The denomina- 
tor of this fraction is determined from Fig. 4. Therefore 
we obtain 
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FIG. 3; Diagram representation of the correlation function Aa-^i . Double solid lines depict the renormalized conduction electron 
propagators Gccr'i'^ ~ t \X)- The arguments of irreducible functions are supposed arranged in the clock wise direction. 
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FIG. 4; Connected vacuum diagrams of the second, fourth and sixth order of perturbation theory. 
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(18) 



The thermodynamical potential becomes equal to 

A 

dX' 



/3 



The expression (19) for thermodynamical potential con- 
tains additional integration over the interaction strength 
A and is awkward because of it. As was proved for non 



^ ^ ^0 + -s /_^ /_^ / -^Gcri'i'^l^ )'Scr[iLu\X ). (19) correlated many-electron system by Luttinger and Ward 



A' 



From this equation we obtain 
dF 2 



P'^l this expression can be transformed into a much more 
convenient formula. 



A—— = — yjzZ G!aiiuj\X)^aiiuj\X). (20) We consider the following expression: 
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FIG. 5: The result of integration over the interacting strengtli of the vacuum diagrams. 



Y = -IY.Y1 cxp(ic^O+){ln[G° (iw)E,(zw|A) - 1] + G,(zc^|A)E,(zc^|A)} + Y\ 



(21) 



which is the generaHzation of the Luttinger-Ward equa- 
tion for strongly correlated systems. Here Y' is the sum 
of closed linked skeleton diagrams with full G„ function 
as a contribution of conduction electron lines. 

On the Fig. 6 are depicted some of simplest skele- 
ton diagrams. These diagrams depend on the interac- 
tion strength A not only through the factors in front of 
each diagram but also through the full Green's function 

Ga{lUj\X). 

From equations (12), (14) and (21) we obtain 
f)Y 1 



dY' 



/? 



dY.„{iu}\Xy 

where, from Fig. 3, 6 and definition (14), it follows that 
dY' A^ S„(iw|A) 



(22) 



^A.(zc.|A) = ^" 



(23) 



dOaiMX) jS-'"-^'-'--' p 

As a result we obtain the stationary property with re- 
spect to changes of the mass operator: 

dY 



9E^(iw|A) 



0. 



(24) 



Now we shall find the quantity ^^ 



-TTT- ■ By the stationary 
property of Y we can ignore the dependence of S^. and G^- 
on A and take into account only the explicit dependence 
of A in y, depicted on the Fig. 6. From this figure it is 
easy to obtain: 



^9A l^- 
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^= «EE^-(HA)S.(*w|A). (25) 
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From equations (20) and (25) we obtain 



OF 
^9A 



dY 



(26) 



The consequence of this equation is the solution 



F{X) = Y{X) + CONST. 



(27) 



For A = we have r(0) = and F{o) = Fq. Therefore 

CONST = Fq. 

The final result has the form 

F = Fo+Y. (28) 



IV. CONCLUSIONS 

The thermodynamic potential of a strongly correlated 
system described by the Anderson impurity model has 
been calculated. We have formulated a new diagram- 
matic technique for fermions with strong correlations and 
determined the correlation function of localized electrons 
and mass operator of conduction electrons. For the con- 
duction electrons this operator coincides with the cor- 
relation function of the impurity electrons. A Dyson's 
type of equation for the one-particle propagators of both 
subsystems, of conduction and impurity electrons, has 
been established. Within our diagrammatic technique 
we first obtained an exact expression for the thermody- 
namic potential as a product of the full propagator G^ of 
the conduction electrons and its mass operator So-, then 
a Luttingcr- Ward-type t^'^l of identity based on the sta- 
tionary property of the potential was established. The 
expression for the thermodynamic potential so obtained 
could be very useful to calculate in a systematic way all 
thermodynamic quantities (e.g. specific heat) of strongly 
correlated electron systems. 
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FIG. 6: Closed linked skeleton diagrams. The double solid lines correspond to full propagators Gcr{ii^\X) of conduction 
electrons. The rectangles correspond to the correlation functions of the correlated electrons. 
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